SOME RESULTS ON INHOMOGENEOUS DISCRIMINANTS 



MARIA ANGELICA CUETO AND ALICIA DICKENSTEIN 

Abstract. We study generalized Horn-Kapranov rational parametri- 
zations of inhomogeneous sparse discriminants from both a theoretical 
and an algorithmic perspective. We show that all these discriminantal 
parametrizations are birational, we prove some results on the corre- 
sponding implicit equations, and we propose a combinatorial algorithm 
to compute their degree in the uniform case of (co)dimension 3. 



1. Introduction 

Given a configuration of n lattice points A in Z^" 1 , let Fa = YlaeA Xa ^ a 
denote the generic polynomial in d— 1 variables (ii, . . . , td-i) with exponents 
in A. Under certain general conditions, Gelfand, Kapranov and Zelevinsky 
|15j showed that there exists an irreducible polynomial with integer coeffi- 
cients Da = Da{x) in the vector of coefficients x = (x a : a € A) (defined 
up to sign), called the ^4-discriminant, which vanishes for each choice of 
coefficients c for which Fa and all its partial derivatives have a common 
root in the torus (C*) rf_1 . The ^4-discriminant is an affine invariant, in the 
sense that any configuration of points affinely isomorphic to A has the same 
discriminant. 

Given A, we form the dx n integer matrix (also called A) whose first row 
consists of ones, and whose columns are given by the points (l,a) for all 
a E A. The kernel of this matrix expresses the affine dependencies among 
the given configuration of points. Let B = (pij) € ^ nx ( n_d ) be a matrix 
whose column vectors are a basis of the integer kernel of the matrix A, i.e. a 
Gale dual of A. B is of full rank and its maximal minors have g.c.d. qb = 1. 
Since the first row of A is the vector (1, . . . , 1) € Z n , the row vectors of B 
must always add up to 0. We say that B is regular when this condition is 
satisfied. 

Set m = n — d. The A-discriminant Da is ^4-homogeneous, i.e it is quasi- 
homogenous relative to the weight defined by any vector in the row span of 
A. Therefore, "taking out these homogeneities" we get a polynomial Ab in 
m variables y± , . . . , y m which is in fact the implicit equation of (the closure 
of) the image Sb of the rational Horn-Kapranov parametrization |15 [ 120 ( 112] 



Both authors are partially supported by UBACYT X042, CONICET PIP 5617 and 
ANPCyT PICT03-20569. 

1 



2 



MARIA ANGELICA CUETO AND ALICIA DICKENSTEIN 



defined by 

n 

(1) yj = JJ(6ii«i H h b im u m ) b ^ , j = 1, . . . ,m. 

i=l 

Extracting the homogeneities means the following. As we said, there exists 
v E Z d such that all monomials c v that occur in Z}^ = ^2, v d v x v satisfy 
A ■ v = v. Therefore, for any v$ such that A ■ vq = v, 

D A {x) = x u ° ^d v x u - u \ 

v 

where d v E Z \ {0} and v — v§ E ker^(j4). Write each v — vq as a linear com- 
bination of the columns . . . , of £?. Then, there is a Laurent polyno- 
mial Ab(u) in m variables such that up to a monomial, A.b{x vW , . . . , x v( ~ m) ) 
equals .D^(x). In particular, A# has the same number of monomials and 
the same coefficients as Da- 

Via the well known Cayley trick, the computation of mixed sparse resul- 
tants may be reduced to the computation of sparse discriminants |15|IH]. The 
computation of sparse discriminants Da (including all sparse resultants) , is 
then equivalent to the implicitization problem for the parametric varieties 
given by These implicitization problems are thus interesting and very 
hard, since they involve the whole of sparse elimination. 

It is easy to see that when the matrix B E Z nxm defining the parametriza- 
tion (^Q) is not of maximal rank m, then Sb is certainly not a hypersurface. 
This is the case in which the corresponding homogenized discriminant Da is 
just the constant polynomial 1. However, the image of the parametrization 
may fail to be a hypersurface even if B is of maximal rank (c.f. Example 12. 4)1 . 
The classification of such defective configurations is a hard combinatorial 
problem [Hj I12| IT5]. On the other side, the condition of non defectiveness 
can be checked algorithmically. 

When trying to compute discriminants, natural reductions as in [H] lead 
to integer matrices B which are still regular and of maximal rank, but whose 
columns do not necessarily generate a saturated lattice ZB, i.e. such that 
the gcd of its maximal minors takes any value gs E Z \ {0}. On the other 
side, sparse discriminants describe the singular locus of ^4-hypergeometric 
systems of PDE's, while the inhomogeneous discriminants A# describe the 
singular locus of classical Horn hyper geometric differential equations |17j . 
In this setting, it is again natural and necessary to consider matrices B 
with arbitrary gcd gs HH|- Another interest on studying this case stems 
in the precise version we give of a result of Kapranov |2U| . which asserts 
that the rational hypersurfaces they define have an interesting geometric 
characterization. 

In this article we focus then on rational parametrizations of the form Q 
for regular non defective integer matrices B E Z nxm of maximal rank m and 
any value of gcd ^ 0. To emphasize the fact that the lattice generated 
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by the columns of the matrix is not necessarily saturated (equivalently, the 
lattice generated by the rows is not necessarily Z m ), we will call our matrix 
C and we will keep the name B in case gs = 1. The irreducible equation of 
the closure of the image Sc, denoted by Ac E Z[yi, ■ ■ ■ , y m ]> will be referred 
to as a generalized inhomogeneous discriminant. 

We prove in Section|2]that, up to multiplication by an element A € (C*) m , 
the parametrizations a la Horn-Kapranov associated to such matrices C cor- 
respond precisely to those rational hypersurfaces Sc for which the logarith- 
mic Gauss map is birational. In fact, we follow line by line Kapranov's 
proof in 20 and we correct slightly his statement in the sense that the 
condition gc = 1 is not necessary. In particular, we show that all such 
parametrizations are birational, thus recovering and explaining the results 
in [H Section 3.3]. 

We then study in Section El the precise relation between the generalized 
inhomogeneous discriminant Ac which gives the equation of Sc and the 
equation A# of Sb, where the columns of B are a Z-basis of the saturated 
of the lattice ZC spanned by the columns of C. We also give other tips to 
simplify the search for the equation Ac- 

Finally, in Section|I]we focus on the case m = 3. The case of codimension 
m = 2 has been studied in detail in We give a combinatorial algorithm 
for the computation of the degree of Sc m case C is uniform, based on the 
intersection formula ()15l) and the algorithm for the computation of local 
multiplicites at the base points given in 9, Prop. 1.5]. This degree could 
be also obtained via Theorem 2.2 in ^2] which should be made explicit to 
give a dehomogenized version of Theorem 1.3 in that paper. It would be 
interesting to develop intersection formulas which do not require to take a 
common denominator in order to work in projective space, to avoid the new 
base points coming from this approach. We also give several examples which 
show the difficulty of computing the local multiplicities in the general case. 

2. Generalized inhomogeneous discriminants and birational 

Gauss maps 

We start by setting some notations and first consequences of the defi- 
nition of the Horn-Kapranov rational parametrizations ipc associated to 
C-matrices and their closed images Sc- We then show that these para- 
metrizations are proper (i.e. deg(V'c) = 1) and that (torus translates of) 
varieties of the form Sc give all rational varieties with birational logarithmic 
Gauss map. 

2.1. The setting. Given a matrix C € Z nxm , n > m, we denote by 
Ci, . . . , C n E Z m the row vectors of C. Each defines a linear form 

(2) l k (ui, . . . ,u m ) := (Ck , (ui, . . . ,u m )). 

Throughout this paper, we will always assume that C is regular (i.e. its 
columns sum up to zero), and has no zero rows. We associate to the matrix 
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C an algebraic variety Sc in C m which is the closure of the image of the 
rational parametrization: 

f C :C m ->C m ( Ul ,...,u m )^( yi ,...,y m ), 

where 

n 

(3) y k = • • • ,« m ) c " Vfc = l,...,m. 

i=l 

Call / = n - =i/ -™«{o, C! ,,:fc=i,.., m } theleagt 

common denominator of all 

the ykS, and write 

(4) y k = ^, k = l,...m. 

Jo 

Since the coordinates of the parametrization are given as a product of linear 
forms whose exponents sum up to 0, all /, have the same degree dc, where 

n 

(5) d c = -^min{0,c ijfc : k = 1, . . . ,m}. 

i=l 

We can think of the mapping ipc as & rational function between projective 
spaces 

(6) Vc : F m_1 — » F m , 

where tpc = (/o : fi '■ ' • • : fm) is defined outside the base point locus 
Z = V(/o, • • • , fm)- We again denote by Sc the projective variety defined 
by the closure of the image of this map. 

The linear forms l k define a hyperplane arrangement in C m and in P m_1 . 
Let T be a flat in this arrangement, i.e. a linear space defined by the 
vanishing of a subset {Z^, ... ,1^} of the given linear forms. Denote by C{F) 
all linear forms lj vanishing on i.e. all linear forms lj which lie in the 
linear span of . . . ,h r }. Note that a given f k vanishes on T if and only 
if it contains a linear factor from C{J-). A given flat T of the arrangement 
will be called basic if all of /o, ... , f m vanish on T. We then have 

Lemma 2.1. The base point locus Z equals the union of all basic flats. 

Remark 2.2. By taking out any common factor from /o, • • • , f m > we can 

always assume that codim(iJ) > 2. Then, in the case of codimension 3 that 
we will study in Section^we can assume without loss of generality that the 
number of base points is finite. Note however that their local structure can be 
very complicated. When m > 3, the base point locus variety has in general 
positive dimension. 

Recall that C is called defective if codim(S'c) > 1. We easily have 

Lemma 2.3. Let C € Z nxm be a regular matrix of rank r < m. Then 
codim(S'c) > m — r, so C is defective. 
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Proof. Assume C\, . . . ,C r are linearly independent and write any other 
= Si=i K.k^k, j = r + 1, . . . ,n. It follows that yi, ■ ■ ■ ,y m can be writ- 
ten as homogeneous rational functions of Zi, . . . , l r of degree 0, and so they 
can be written as rational functions of the r — 1 variables h/l r , ■ ■ ■ , l r -i/l r . 
Therefore, the codimension of Sc is at least m — (r — 1) = m — r + 1, as 
claimed. □ 

Example 2.4. As we have remarked in the Introduction, the converse to 
Lemma 12.31 is not true and the classification of defective matrices involves 
subtle combinatorial questions. As a simple example of a full rank defective 
matrix, let n = 2n' be even and m < n' arbitrary. Pick any set of integer 
vectors Ci, . . . C n > whose Q-linear span is Q m and take C n i + k = — Cfc f° r all 
k = 1, . . . ,n'. Then, it is straightforward to check that tpc is a constant 
map. 

We will assume from now on that C has full rank m, so that the gcd of 
its maximal minors gc is non zero. 

2.2. Birationality of ipc- Given an algebraic group G with Lie algebra 
Q, let l g : h i— > gh denote the left translation by an element g G G. Let 
S C G be an irreducible algebraic hypersurface. The left Gauss map of 
the hypersurface S is the rational map 75 : S — > P(£7*), taking a smooth 
point y G S to the hyperplane dQ^XTyS) C T e (G) = 0. The case which 
we are interested in is G = (C*) m . So Q = C m , and left translation is the 
usual coordinatewise multiplication map by an element in the torus. If S is 
a hypersurface defined by a minimal equation (A = 0), then the Gauss map 
is simply 

(7) i{y) = (ml— (y) ■ ■ ■ ■ Vm^ L {y)) , 

oyi oy m 

mapping a regular point y £ S to a projective point in P" 1 ^ 1 . Assume 
S* = S n (C*) m is non empty. Geometrically, the map 75 corresponds to 
looking at the image of S* via the map log(y) = (log(yi), . . . ,log(y m )) and 
then considering the Gauss map of log^*). 

The following theorem is essentially due to Kapranov [2111 • Our con- 
tribution lies in making the statement precise by removing the incorrect 
hypothesis about the gcd of the matrix defining the parametrization. We 
will show in Section 01 why a quick first thought may lead to conjecture that 
this gcd should equal 1 (cf. E,emark l3.2j) . 

Theorem 2.5. Let S C C m be an algebraic irreducible hypersurface. 

The Gauss map 75 : S — ■* P m_1 is birational if and only if there exist 
a non defective and regular integer matrix C € Z nxm of full rank, and a 
constant A G (C*) m such that S = A • Sc, i-c. S is a torus translate by A of 
a generalized inhomogeneous discriminant hypersurface. 

Moreover, in this case, A • ipc ^ birational and the logarithmic Gauss map 
7s is its inverse. 
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As we said in the Introduction, the proof is exactly Kapranov's original 
proof, except that one can check that the hypothesis gc = 1 is not used 
in the "if" direction and that in the "only if" direction his last statement 
about the fact that gc needs to be equal to 1 for ipc to be of degree 1, is not 
true. For the convenience of the reader, we sketch only the "if" direction, 
i.e. the proof that for any C as in the statement, the logarithmic Gauss map 
75 c is indeed a birational inverse to ipc- 

We first give a simple example. 

Example 2.6. [The discriminant of a generic univariate cubic polynomial] 
The vl-discriminant associated to the 2x4 matrix 



A 



1111 
12 3 



B 



i.e. the discriminant Da(xq, xi,X2,xs) of the generic polynomial Fa = xq + 
x\t + X2t 2 + x^t 3 equals 

Da{x) = -21x\x\ + 18x 3 x x 2 xi + x\x\ - Ax\x§ - Ax^x\. 

Equivalently, {Da = 0) is the dual variety of the projectively embedded 
toric variety parametrized by monomials with exponents in A, i.e., the well 
known twisted cubic. 

Consider the following choice of matrix B: 

( 1 2\ 

-2 -3 
1 

Note that gs = 1, so that the columns of B are a basis of the integer kernel 
of A. Calling Iq(u,v) := u + 2v , l\{u, v) := — 2u — 3v , l2(u,v) := u and 
ls(u,v) := v, the parametrization ipB equals 

yi ■■= ^ 

llh 

Its closed image is the hypersurface Sb = = 0}, where 

A B (yi,y 2 ) = -4y 2 - 27y| + y\ + 18y 2 yi - 4yf, 

which can be computed in Maple as 

—normal (re suit ant (if * yi — 1$ * if * V2 ~ Iq * hi u ) I 

(the factor v 6 appears because Maple computes an affine resultant instead 
of an homogeneous one with respect to (u : v)). Equivalently, it could be 
computed as the dehomogenization 

As (2/1,2/2) = D A (1, 1,2/1,2/2). 
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Conversely, up to a monomial Da equals Ab(xqX2/x 2 , XqX^/x 3 ). The asso- 
ciated logarithmic Gauss map equals 7_b = (71,72), with 

71 := I2y\ - 2y\ - 18 yi y 2 , 72 := -4y 2 - ^2 + 18j/ilte- 

If we compute (21 (u, u), .z 2 («, «)) := (7_b iPb)^, v), we obtain: 

2l = A(u + 2v) 2 u(u 3 + 9u 2 v + 27uv 2 + 27v 3 )/(2u + 3vf , 

z 2 = 4(n + 2v) 2 u(u 3 + 9u 2 t; + 27nu 2 + 27t; 3 )/(2n + 3?;) 6 . 



Note that we do not recover (u,v) but it only holds that z\/zz = u/v, or 
(zi : z 2 ) = (u : v). 

Consider now the following matrix C: 



C 



Note that A ■ C = but that g c = 3. In fact, C = B ■ M, where M is the 
square matrix: 

-3 
2 1 



/ 


1 


2 \ 







-3 




-3 





V 


2 


1 / 



M 



Calling Lq(u, v) := u+2v , L\{u, v) := —3v , Lz(u, v) := —3u and L%(u, v) 
2u + v , the parametrization equals 



2/1 •— - pr- 

ft =- ^ 

Its closed image is the hypersurface Sc = {Ac = 0}, where 

A c (yi,y 2 ) = -12962/ 2 y 3 - 874% 2 V - 19683y|y 3 + y 2 yi + 4698yiy 2 
-64y 3 - 64y 3 + 24y 2 V + 24y 2 y 2 - 1296y 3 yi - 8748y 3 y 2 . 

The corresponding logarithmic Gauss map jc(y) composed with ipc(u,v) 
gives back (r(u,v) • u,r(u,v) • v), where r denotes the rational function 
r(u,v) = -4/729(w + 2vf{2u + u) 3 (16v 4 + 84uu 3 + U3u 2 v 2 + 84u 3 v + 
16ti 4 )(i> +u) 3 / (u 6 v 6 ). We will study in Section|I]the precise relation between 
the irreducible polynomials Ab and Ac- 

Proof of the "if" part in Theorem \2.5\ Let C be a regular non defective n x 
m integer matrix, a point A £ (C*) m in the torus, and consider the map 
ip' c := Xipc- We need to show that the logarithmic Gauss map is its bira- 
tional inverse. Denote by A an irreducible equation of its closed image. The 
principal observation is that the Jacobian matrix of log(^) is symmetric 
since 

a n 
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Moveover, a straightforward computation shows that for any point u in the 
preimage of the torus, the Jacobian matrices J(ip'c) and J(log(^' c ) have the 
same rank since J(\og{ijj' c ) = J(iP'q) ■ D, where D is the diagonal matrix 
with diagonal entries the multiplicative inverses of the coordinates of ip' c . 
This rank is equal to m — 1 by our hyphotesis that C is non defective. Now, 
on one side, implicit partial differentiation of the equality A(ip' c (u)) = 
implies that the vector 7c(y) lies in the kernel of the transposed Jacobian 
matrix J(log(ip' c )) t for any y in the image of ip' c . On the other side, since 
the coordinates of ip' c are homogeneous forms of degree 0, it follows from 
Euler's formula applied to the coordinates of log(ip' c )) that any point u in 
the preimage of the torus lies in the kernel of J(\og(ip' c )(u). Then, u is 
proportional to 7c(V , c( M ))> when this vector is non zero. □ 



3. Monomial changes of coordinates and factorizations 

Given any matrix C G ^ nxm with gc ^ 0, we proved that the parametri- 
zation ipc is birational if C is non defective (i.e. when Sc has codimension 
1). Indeed, the converse is trivially true. The defectiveness condition can 
be checked by computing the generic rank of the Jacobian matrix J(ipc){ u )- 
The (m — 1) x (m — 1) minors of this matrix are rational functions, and 
with probability 1 at least one of them will be non zero at a generic point 
u when C is non defective. The tropical approach in provides another 
algorithm to check this property. 

It is obvious that given a regular non defective matrix C G Z nxm , we can 
replace all row vectors in C lying in the same one-dimensional flat T by their 
sum, without essentially changing the coordinates of the parametrization tpc 
except for constants (if the sum gives the zero vector, we keep the constants 
but we don't keep a zero row). So, we can always start simplifying the 
problem by making this replacement. Note however that we could change 
the gcd, so this operation could lead from a Gale dual matrix with gcd equal 
to 1 to a matrix with arbitrary (non zero) gcd [Sj. 

Assume C and B are non-defective n x m regular integer matrices of 
full rank m such that the lattice TLB generated by the columns of B is the 
saturated of the lattice ZC. Then, there exists a square matrix M G ^ mxm 
with determinant equal to ±gc such that C = B-M. as in Example l2.fil The 
lattice ideal /(ZC) (in n variables) is radical with primary components, 
which correspond to torus translates of the toric variety defined by the lattice 
ideal I {TLB) We will see in Theorem 13.41 how this is reflected in the 

precise relation between the irreducible m-variate polynomials A# and Ac- 

We begin by recalling Lemma 6.3.1 in [Jj, which shows the relation be- 
tween the Horn-Kapranov parametrizations associated to two regular integer 
n x m matrices C\, C2 with C\ = C2 ■ M for a square matrix M. 
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We associate to M = (My) € Z mxm a linear map A M : F"™" 1 P" 1-1 : 
(8) A M {u) = Mi^Uj : . . . : M ™,jUj) = M ■ u l , 

3=1 3=1 

and, denoting by M^ l \ . . . , the column vectors of M, the (multiplica- 

tive) monomial map om ■ (C*) m — » (C*) m : 



m m 

(9) a A /(y) = (llvi ',■■■,11% ) = (y >••••!/ 

i=l i=l 



For any given m x m matrices Mi,M2, it clearly holds that A^.jv^ = 
Ami o A M2 and 

(10) «A/i-Af 2 = a M 2 &Mi ■ 

Lemma 3.1. The following diagram 

1 A M 1 
jpm-1 >. jpm— 1 

I I 
V'Ci l l i>c 2 

ri m ^ (cr 

is commutative. 

The proof of Lemma l3,ll is a straightforward verification. 

Remark 3.2. Assume C2 has gcd equal to 1, and call C = C\,B = Ci- 
Then, | det(M)j = gc- Suppose that we didn't know Theorem \2.5\ but instead 
we suspected (or proved) that i/jb is birational. From the equality ipc = 
olm otpB Ajtf, where Am is birational and olm is a gc to 1 mapping, one is 
tempted to deduce that t/jc is also a gc to 1 mapping. But indeed, we have 
already proved that it is birational. The explanation is given in the next 
lemma. 

Lemma 3.3. With definitions and notations as in and 0), the restric- 
tion of the map olm to the zero set of Ac 2 defines by corestriction a birational 
map 

olm = «Af| (Ac , 2=0) : (Ac 2 =0)—* (A Cl = 0), 

where as before C\ = C2.M and C\,Ci are non defective regular integer 
matrices. 

Proof. The image of ipc 2 is dense in (Aq 2 ) = 0, and so the image of ipc 2 °Aj\/ 
is also dense. For any point y of the form y = ipc 2 (A-M(u)), its image by 
ajf equals the point ipCi{v>)- Therefore, it lies in (Ad = 0). Thus, we have 
a rational map olm '■ (Ac 2 =0) ™ » (Ad = 0) which has to be 1-1 by 
Theorem 12.51 applied to both ipc 2 and ipd- ^ 
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We now present the relation between Ad and Ad- Note that since 
the corresponding Horn-Kapranov parametrizations are given by rational 
forms with rational coefficients, we can assume that these polynomials have 
integer coefficients and content 1. They are thus defined up to sign, but we 
will usually omit this sign in the notation. 

Given an integer square matrix M of size m, denote by Gm the multi- 
plicative group 

(11) G M := {e E (C*) m = a M (e) = (!,..., 1)} 
with the induced coordinatewise multiplication. We then have: 

Theorem 3.4. Let C%, C2 are non defective n x m regular integer matrices 
such that C\ = C2.M. There exists v in the lattice ZM spanned by the 
columns of M (or equivalently, verifying e v = 1 for all e E Gm ) such that 

(12) A Cl oa M (y) = y v II A c 2 (e-y). 

e£G M 

Before giving the proof, we revisit Example 12.61 

Example 3.5. [Example 12 . 61 cont .] With the notations of Example 12.61 the 
group Gm consists of those (£1,62) E (C*) 2 such that e\ = 1, £2 = 1 and so 
we have 

A C (4^2) = 4 IT ^B(ey u y 2 ). 

Of course, we can move the factor yf to the left hand side to produce an 
equality in the polynomial ring C [2/1,2/2] but this is really an equality over 
the Laurent polynomial ring C[2/i,2/2 1 ]' The exponent vector v = (—9,3) 
equals 3 times the difference of the columns of M . 

Proof. Given any point y E (C*) m in the image of ipCa (thus, in a dense 
subset of (Ad = 0)), using that Am is an isomorphism and Lemma l3~Tl we 
can write y = ipc 2 (A-M(u)) for some u E p m_1 and so € (Ad = 0). 

For any e £ Gm we have that om(2/) = &m( £ ' y) an d therefore it also 
holds that Ad(e • y) = 0. Reciprocally, pick any point y G (C*) m such 
that dM(y) lies in the image of ipCi- Then, there exists u E p m -i SU ch 
that «m(V ; C2 (A-m(u))) = tpd(u) = a,M{y)- Therefore, there exists e E Gm 
such that y = e~ 1 tpc2(A-M(u)), and then Ac 2 (£ ■ y) = 0. By density and 
properness arguments, we deduce that that 

(Ad o aM (y) = 0) n (c*r = U ( A ^( £ • y) = °) n ( c *) m 

eeGm 

Observe now that, as a consequence of Lemma 13. 3| the irreducible po- 
lynomials Ad(e • y) with e varying in Gm, are pairwise coprime. In fact, 
if Ad(e • y) is proportional to Ac 2 (e' • y), with e,e' E Gm, then writing 
S = e' ■ e -1 , we have that a point y E (Ad = 0) n (C*) m if and only if the 
point 5 ■ y E (Ad = 0) n (C*) m . Since cum is birational, we deduce that 
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5 is the unit element in Gm, i-e. that e = e' . Then, we deduce from the 
Nullstellensatz that there exist positive integers n e , and a unit qy v (where 
q is a constant and v £ Z m ) in the Laurent polynomial ring Cfy^ 1 , . . . , y^ 1 ] 
such that 

(13) A Cl (a M (2/)) = qy v ]J A C2 {e-y) n *. 

eGGm 

Substituting y i — ^ <5 - y for any 5 € Gm m the above factorization we get 
A Cl («A/(y)) = A Cl (a M (S-y)) = q5 v y v ]J Ac^e-y)"-'- 1 . 

By uniqueness of the irreducible factorization, it follows that all n £ are equal 
to some JVeN and that moreover (5 11 = (1, . . . , 1) for all 5 € Gjy. 

It is clear that this last property holds whenever v lies in ZM. To prove 
the converse, assume that 8° = 1 for all 5 € Gju- To see that v £ ZM we 
write M in its Smith normal form: 



(14) M = U ■ 



( di 
\ d m 



•V, 



=D 

where U, V 6 r E mxm are invertible over Z and di | cfe I • • • | dm in By 
the composition formula (|1U|) . it is enough to prove a similar result for U, D 
and V. This is obvious for U and V since ZU = ZV = Z m and it can be 
easily proved for the diagonal matrix D by choosing elements in Go of the 
form (1, . . . , Wk, . . . , 1) with Wk a primitive root of unity of order dk, for all 
k = 1, . . . ,m. 

Assume that all n £ = N > 1 and differentiate both sides of equation ()13|) . 
Since the Jacobian of au is invertible at any point of the torus, we deduce 
from the Chain rule and the fact that olm is birational, that the Jacobian 
of Acj vanishes along (A^ = 0). But this contradicts the irreducibility of 
Ac 1; so N must equal 1. 

Finally, we show that q = dbl. As A^ has content 1, we need to show 
that the product 

P(y):= [] A C2 (e-y) 

eeG M 

has integer coefficients and content 1 too. First, note that if we write M 
in its Smith Normal Form (|14)L it is enough by formula (|10j) to prove the 
result for the factors U,V,D. This is obvious for U and V. We can further 
decompose D as a product of diagonal matrices, all of whose diagonal entries 
are equal to 1 except for a single entry which is a prime p. So, assume that 
D = (dij) is the diagonal matrix with da = 1 for all i = 1, . . . , m — 1, and 
dmm = P- But then, the coefficients of P are symmetric polynomials with 
coefficients in Z[sx, . . . , s p ], where the elementary symmetric 

functions on the p-th roots of unity. Since all equal either 1,0 or —1, we 
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deduce that P 6 Z[yi, . . . , y m -i][y m ]- In fact, there exist a polynomial Q 
with the same coefficients, such that P(yi, . . . , y m ) = Q(yi, • • • , Vm-l, Vm)- 
Moreover, for every fixed values of the first m — 1 coordinates, the roots 
of Q are the p-th powers of the roots of Ac 2 in the last variable y m , so 
that one can trace recursively the relation between the coefficients of Q and 
the coefficients of Ac 2 (which has content 1) to deduce that the gcd of the 
coefficients of Q is also equal to 1. 

A different argument to prove that q = ±1 is the following. Assume again 
that we have a diagonal matrix with all diagonal entries equal to 1, except 
for a single entry which equals a prime number p. Fix a primitive p-th root 
of unity w and take any monomial degree ordering -<. Call 6 7 y 7 the leading 
term of Ac 2 (y). Then, the leading term of P is tiyW^ m ^t=i 8 )y#T = ±\ b^y 9 ^. 
Since Ad £ Z[y], then q G Q. But as all the coefficients of P lie in Z(w), 
and q G Q we deduce that P € Z[y]. In fact, g = 1/s where s =cont(P). 
Assume s / il and let a be a prime dividing s. Suppose first that a / j). 
Since a divides the content of P(y) G Z[y] we have that P = in the 
extension field 7j a {w). Therefore, one of the factors of P must be zero. But 
given that the content of Ac 2 is ±1, this cannot happen. If a = p, first 
reduce Ac 2 mod p and then look at its leading coefficient, which we call 
6 7 y 7 . Then, the coefficient of the monomial y 91 in P is not divisible by p, a 
contradiction. □ 



The moral of the factorization provided by Theorem \'AA\ is that when 
trying to compute a generalized discriminant Ac it is possible to compute a 
"simpler" discriminant A# where ZC C "LB and then reconstruct Ac from 
(|12ft . Even if the monomial y v is in principle unknown, its function is to 
clear denominators without introducing extra monomial factors. The group 
Gm can be computed via the Smith Normal Form l)14|) of M. One easier 
way to recover Ad from Ac 2 is the following. Substitute y = axd](M)i z ) 
in lfT2|) and denote g := gc 2 - The key point is that au a Adj(M)( z ) = z9 ■ 
Furthermore, since we know that the exponent v lies on the lattice spanned 
by the columns of M, we deduce that the specialized product on the left 
hand side must be a Laurent polynomial in the variables zf,...,Zm- So, we 
only need to divide all exponents in the resulting polynomial by g to recover 

Ad(*)- 

In fact, this is also useful when dealing with saturated lattices ZP, i.e. 
when trying to compute ^-discriminants. Given a regular n x m non defec- 
tive matrix B, we can first look for a reduced basis of ZP using the LLL- 
algorithm |18j . available at any Computer Algebra System, and then put 
these reduced generators as the columns of a matrix B' . Write B = B' ■ M 
with det(M) = 1. Then Gm consists of the single element (1, . . . , 1) and so 
(|T2|) reduces to the equality 



A B {y) = y v A B ,{a M -x{y)) , 
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which can in fact be easily proved since the homogenizations of both dis- 
criminants give the same discriminant Da (Here, B and B' are Gale duals of 
A). Since the coefficients in the parametrization i^b 1 are smaller, the implicit 
equation A b> can be obtained using standard elimination techniques in cases 
in which the systems would crash when trying to compute A# directly. 
Here is a simple example. 

Example 3.6. [Example l2 . 6l cont .] Recall that for the matrix B we proposed 
in Example 12. 61 the inhomogeneous discriminant had degree 3. Consider now 
another choice B' of a Gale dual of A whose entries are integers with larger 
absolute value: 



/ 


-5 


" 3 \ 




13 


8 




-11 


-7 


V 


3 


2 / 



If we try to compute the inhomogeneous discriminant A b 1 using the resul- 
tant formula, we get a common factor in the coefficients equal to 3 88 . After 
dividing by this quantity, we recover the following polynomial of degree 16: 

A B ,(y) = -27yf + lSyfy? - ±y\y\ - 4y 2 3 y 8 + yf. 

Note that in general, the degree of a dehomogenization of a sparse dis- 
criminant Da can be as large as wanted. 

4. The degree of Ac and the computation of local 

multiplicities 

Let C E ^ nxm b 6) as before, a regular non defective integer matrix with 
no zero rows. We set m = 3 and we assume w.l.o.g. that the variety of 
base points is finite. Base points p? £ Z are indexed by basic flats T as in 
Section [2j In this section, we concentrate on the algorithmic computation of 
the degree of the generalized homogeneous variety Sc, based on the following 
well known intersection theory formula ^1] 

(15) d 2 c = deg(xp c ) deg(S c ) + ^ ejr, 

T basic 

where denotes the Hilbert-Samuel multiplicity of p? [TTH l2"l IHj. 

Since we have an easy formula for dc in © and we know that deg{ipc) = 1 
by Theorem 12 .5| we would need to compute the Hilbert-Samuel multiplici- 
ties. In fact, this is a delicate notion and there is no efficient deterministic 
algorithm for the general case (cf. |16| 121 j for a Grobner /Standard bases 
approach). We will start by recalling the definition of multiplicity, together 
with some known properties. In particular, there is a probabilistic algo- 
rithm, which reduces the problem to the tractable case of local complete 
intersection. We will present however several examples, which show that 
in general one cannot expect that the base points that occur in discrimi- 
nant parametrizations are local complete intersections, or even almost local 
complete intersections with one more generator. It follows that the known 
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algorithms to find the implicit equation do not work in principle in these 
cases p] H] . One way out would be to compute the Newton polytope of Ac 
using the results in ^2j, and then compute its coefficients via an efficient 
interpolation. 

We show in Proposition ^. 31 that when C is uniform, i.e. when all its max- 
imal minors are non zero, one can use the combinatorial algorithm from [2] 
to compute the Hilbert-Samuel multiplicities. We then turn things upside- 
down in Corollary 14.81 to compute the dimension of the local vector space at 
the origin of d sparse polynomials in d variables with generic coefficients. 

There are several ways to define the algebraic multiplicity of a base point. 
Our definition follows [T^j . We will state the results for the case of dimension 
two, but they hold with the obvious changes in any dimension. We refer to [5] 
for more details on the definitions and examples in this section. 

Let p = p?r g Z. Consider the Noetherian local ring A p := Op2 p and the 
localized base point locus ideal I p := (fo, fx, fa, fs)A p . The Samuel function 
of A p with respecto to I p is defined as: 

X I X p (r) = l(A p /i;+ 1 ) forallrGN, 

were Z(_) is the length function of I p as an Ap-module, that is, the length 
of a composition series of the module. Since we are working over the alge- 
braically closed field C, this length coincides with the vector space dimen- 
sion dime (Ap/ I p +1 ). The Samuel function is polynomial for large values 

of r, that is, there is a polynomial PS^(X) in Q[X] (which takes integer 

values over Z) such that we have PS^(r) = XA p ( r ) f° r r >> 0. Moreover, 
this polynomial has degree 2 and its leading coefficient is e/2! with e £ No- 
Then, the local multiplicity ejF of the base point locus at p?r is defined as 
ejr ■= e, i.e. 2! times the leading coefficient of the polynomial PS^. 
Therefore, 

dimcCV^ 1 ) • 21 
ejr = lim ^ . 

r^oo f 

When the base point p = p^r is a local complete intersection, i.e., when the 
ideal I p admits two generators, then ejr is just the vector space dimension 
Ap/Ip of the local quotient. This dimension can thus be computed algo- 
rithmically via a standard basis computation using a local order -<, and 
counting the number of monomials not in in~>(I p ). Even if this is not the 
case, we always have the following probabilistic approach to compute the 
local multiplicity. 

Consider the ideal J p generated by 2 generic linear combinations of the 4 
generators: 

J P :=< v° fo + v°xfx + v° 2 f 2 + v° 3 h , vlfo + v\h + v\h + ^3/3 > , 

with vj £ C. Then J p is generically a complete intersection inside I p (and 
a reduction ideal of I p ). Thus, with probability 1, we can compute = 
dim c (A p /J p ). 
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As a corollary, we always have the inequality > dime (A p /J p ), so that 
in any case 

deg{A c ) < d 2 c - ^ dim c (A pr / I PT ) . 

T basic 

On the other side when / is a monomial ideal, there exists a combinatorial 
way of computing this multiplicity, as stated in [Hj . We reproduce this result 
below. 

If p = (1 : : 0) is a base point (which we can assume after a translation) 
and the localized ideal I p is monomial, we have the following algorithm to 
compute the Hilbert-Samuel multiplicity e p at p: 

Algorithm 4.1. Computation of Hilbert-Samuel Multiplicities for the mono- 
mial case and m = 3. 

• Set xq = 1 and let I p be the specialization of the ideal I p . 

• Compute the convex hull C of the exponents of the bivariate mono- 
mials in I p . 

• Then: e p = 2! • FoZ(Nq n C) equals the normalized volume of the 
complement K, of C in the first orthant. 

Note that in the very simple case in which the ideal is both monomial 
and a complete intersection, generated by {x™ 1 , x™ 2 }, the local multiplicity 
mi x m2 equals both the normalized volume of the triangle K, with ver- 
tices (0,0), (mi, 0), (0,7712), which is the complement in Nq of the convex 
hull of the staircase of the ideal, and the number of standard monomials 
{x^x^ 1 < ki < mi,i = 1,2}, which is the dimension of the quotient by 
the ideal. 

Remark 4.2. The Volume Algorithm ^. 1\ does not work for general ideals. 
If fact, it might seem reasonable to expect the same algebraic multiplicity for 
Ip and for any initial (monomial) ideal in~>(I p ) with respect to a local order, 
but this is not in general the case. We illustrate this issue in Example \4-5\ 



Proposition 4.3. Assume C S Z nx3 is uniform. Then C is non defective 
and each local ideal I p becomes a monomial ideal modulo a linear change of 
coordinates. So, the degree of the generalized discriminant surface Sc can 
be combinatorially computed using formula \15\) and Algorithm ^. 1\ 

Proof. The fact that C is non defective follows from > an d more precisely 
from [SJ Section 5]. Since any flag is generated by only two linear forms /j, lj, 
after dehomogenizing and localizing, the four generators of each local ring I p 
are products of powers of Zj, lj, i.e. they are monomials in two independent 
variables Zj , lj . □ 
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Example 4.4. We first give a very simple example to illustrate Proposi- 
tion E21 Let C be the uniform matrix: 



C 



( 


2 


1 


M 




-2 


-1 


-2 




1 


1 





\ 


-1 


-1 


- 1 / 



We read from the parametrization tpc that 
fx = l\h , h = hhh , /3 



/V 



ip4 



/o 

1) and p' 



P{1,2} 



There are two base points: p = P{i^} = (—2 : 1 

(1 : —2 : 0). The localized ideals are the following monomial ideals I p 
(ll,h,lf,l4) = (htU) and I p / = (Zf , lifo, if, Zf) = (ifjihJl) (in variables 
Zi , Z4 and l±,l2 and not in the u variables). The first ideal is moreover a 
complete intersection. It is straightforward in this case to compute the 
multiplicities: e p = 1 and e p / = 4, while the dimension of the local ring at p' 
equals 3. Thus, deg(Sc) = 3 2 — 1 — 4 =4. Indeed, in this case gc = 1, so 
that C is a Gale dual of the matrix A £ Z lx4 with all four entries equal to 1. 
So, the homogeneous ^-discriminant Da{x) = x\ + X2 + £3 + £4. Then, we 

can obtain A# by dehomogenizing Da- From the equations y\ 



X\X-j 
X2X4 ' 



V3 



x\ 
— 



xixo, 



read from the columns of C, we get that 



D A (x) 



X\ 



1+ X2 + X3 



X\ 
X\ 



X\ 



1 + y_i + y_i 

yi y-i 



+ 



y\ 



Clearing denominators, we get the equation 

Ac(y) 



y\yz + yw2 + y\ + y\y\, 



of degree 4, as predicted. In fact, this same procedure holds for any matrix 
B of size (m + 1) X m and gs = 1- The associated discriminant A# has 
(m + 1) monomials and all coefficients are equal to 1. 

We address now a more complicated example, which will help us illustrate 
several features. 



Example 4.5. 



Consider the matrix 
/ 



C 



V 



\ 
1 




-2 

1 J 



Observe that the first and the third rows of C are identical. We have: 
/o = li^h ) /1 = ^5^6 j /2 = ifl^h^Q j /3 := ^ihUl^e- There are seven base 
points: £>{i, 3 ,4} = (0:0:1), P{i |2 ,3,5} = (1:1:0), p {lfi} = (1:1:1), 
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y 



(0,3) 




Figure 1 . Region corresponding to the ideal I p 




C 



u 

Figure 2. Region corresponding to the ideal in^{I p ) 

P{2,A} = ("2 : -1 : 1) , p {2 ,6} := (1 : 2 : 1) , p {4M = (-4 : -2 : 1), and 
P{5,6} = (1:3:1). 

Let's focus on p = Pu^m- The local ideal equals I p = (if ; i| ; ifh), or 
changing the name of the variables, I p = {xf , x|, x\x2)- By the volume 
formula in Algorithm 14.11 we get e p = 10 (see Figure The dimension of 
the local quotient by I p equals 8 < 10. 

On the other side, if we write the linear forms li in the affine coordinates 
(ui,ii2), we look at the generators of I p in the polynomial ring <C[ui,U2[- If 
we consider the local order -<= ds (with U2 ~< u\) in Singular we 
get the following initial ideal: 

in^(I p ) = (u\,u\u 2 ,uiu%uf), 

which, by the same volume algorithm as above, has multiplicity 11 > 10. 
See Figure |2j 

The implicit equation Ac(y) can in this case be easily computed with 
Singular by Grobner bases methods and we get (up to sign) the following 
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polynomial of degree 13: 

-8yfyly 2 3 + Syhhs + 1%| + W00yfy 2 y 2 + 3y iy2 y3 + y\ 
yh3 + 2/3 + 3125yf y\y\ + 27y 2 y 2 + 16yfyf yf - 225y 2 y 2 y 3 - 
-225y 3 y 2 y 3 + 500 yi 2 y 2 y2 + 500y 4 y 3 y 2 + l 6 0yf + 80 

vivivl - ^vlvhl - ^vlvhl + l ^y\yhl - ^ymyj- 

8yl + 80yf yiy 2 • 

Example 4.6. One could try to reduce the computation of multiplicities to 
the case of monomial ideals in the following way. Consider for instance the 
matrix 



C 



( 



\ 



1 

-1 






and the base point p 
we have that 2/3 = — 



= P{2,3,4}- Calling xi = i 4 (l, « 2 , it 2 ), x 2 = k(l,ui,u 2 ) 
(U + ^5), so the local ideal I p is generated by 

Ip = {(xi + x 2 )xix% (x 1 + x 2 ) 3 , x\x 2 ). 

This is not a monomial ideal in these coordinates, but after the linear change 
x\ = u + 3v, x 2 = —u + v, it becomes a monomial ideal with the following 
generators obtained after easy algebraic manipulations: (u 4 v , n 6 , w 3 }, and we 
can compute its multiplicity e p = 18 by means of Algorithm 14,11 However, 
this case is very special and there doesn't seem to be a general pattern about 
when such a change of coordinates is possible. 

For instance, it is possible to prove that the local ideal at the base point 
P{i,4,5} °f the parametrization tpc associated to the matrix 



C 



( 



\ 



3 \ 
2 
1 


-6/ 



cannot be transformed into a monomial ideal by a linear change of coordi- 
nates. 

Example 4.7. Consider the matrix C: 





( 


1 


-7 


-6\ 






-1 


4 


3 


c = 




1 





4 









1 


-1 




\ 


-1 


2 





and the base point p = P{i >2 } = (—1 : —1 : 1). 
the local ideal at p equals I p = (x 8 , y 5 , xy 2 , x 7 y) 
complete intersection. 



Calling x = li, y = l 2 , 
which is not an almost 
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e = dim c 



We end with an application to unmixed sparse polynomial systems. Fix 
an exponent set A = {ax, . . . ,a r } C Ng, with r > d, and consider d generic 
polynomials Fx, . . . , F d with exponents in A and coefficients in C: 

r 

F i( x ) = ! i = l,...,d ; x = (xx, ■ ■ ■ ,x d ) 

with (cj) ., • • • , ( c j)j ^ C- r generic. Then, by Bernstein's theorem, the total 

number of common roots in the torus (C*) d equals the normalized volume 
of A. Assume moreover than on = Aj e« for all i = l,...,d, i.e. that a 
monomial which is a pure positive power of each of the variable occurs 
in the polynomials Fx,... , F d . Using the previous results, we are able to 
compute geometrically their multiplicity at the origin (see also Chapter 
5, § 2.E] for a general version of this result). 

Corollary 4.8. Let A = {ax, ■ ■ ■ , a r } C Nq such that on = i = 1, . . . , d, 

where A, € N. Given generic sparse polynomials Fx, . . . , Fd with exponents 
in A, their multiplicity at the origin 

( C\xx, ...,x d } \ 
\(Fx,...,F d )J 

coincides with the normalized volume of the complement Ka in the first 
orthant of the convex hull of {ax + M> d } U . . . U {a r + M> }. 

Proof. By our hypothesis about A, it follows that the monomial ideal I 
generated by {x ai , . . . , x° r } is supported at the origin G C d . We identify it 
with its localization Iq = I <C[xx, ■ ■ ■ , Xd]o- Note that the localized ideal Jo := 
(Fx, . . . ,F d )o is a generic complete intersection inside the zero dimensional 
ideal Iq . By the probabilistic algorithm (in dimension d) for the computation 
of the local multiplicity e of Jo, we know that e equals the local multiplicity 
of the reduction ideal Jo of Io, and then e = cq. Finally, e coincides with 
the normalized volume of JCa by the corresponding version of Algorithm 14. II 
in dimension d. □ 

Acknowledgments. We thank the Institute for Mathematics and its Ap- 
plications (IMA) in Minneapolis, USA, where this work was completed. We 
are also grateful to Eduardo Cattani, Teresa Krick and Rafael Villarreal for 
very useful discussions. 

References 

[1] L. Buse and M. Chardin. Implicitizing rational hypersurfaces using approximation 

complexes. J. Symbolic Computation 40, 1150-1168 (2005). 
[2] W. Bruns and J. Herzog. Cohen-Macaulay rings. Cambridge studies in advanced 

mathematics, 39. Cambridge University Press, 1993. 
[3] E. Cattani, A. Dickenstein, and B. Sturmfels. Rational hyper geometric functions. 

Compositio Mathematica 128, 217-240, 2001. 
[4] L. Buse, D. Cox, and C. D'Andrea.ImpL'citization of surfaces in P 3 in the presence 

of base points. J. Algebra Appl. 2 (2003) 189-214. 



20 



MARIA ANGELICA CUETO AND ALICIA DICKENSTEIN 



[5] D. Cox. What is the multiplicity of a base point?. Slides from the ex- 
pository talk given at XIV CLA, La Falda, Argentina, 2001. Available at 
http://www.amherst.edu/ "dacox/. 

[6] M. A. Cueto. Tesis de Licenciatura: Multiplicidades y variedades discriminantales. 
Universidad de Buenos Aires, 2006: http://cms.dm.uba.ar/lic/tesis/cueto.pdf 

[7] R. Curran. Toric ideals and discriminants in codimensions greater than two. Doc- 
toral dissertation. University of Massachusetts, Amherst, 2005. 

[8] R. Curran and E. Cattani. Restriction of A-Discriminants and Dual Defect Toric 
Varieties. To appear: JSC, 2006. 

[9] D. Delfino, A. Taylor, W. V. Vasconcelos, N. Weininger and R.H. Villarreal. Mono- 
mials ideals and the computation of multiplicities, in: Commutative Ring Theory 
and Applications (M. Fontana, S.-E. Kabbaj and S. Wiegandl, Eds.), Lectures 
Notes in Pure and Applied Math. 231, Marcel Dekker, New York, 2002, 87-107. 
[10] A. Dickenstein, L. Matusevich, and T. Sadykov. Bivariate hypergeometric D- 

modules. Advances in Mathematics, Vol. 196, Issue 1, 2005, 78-123. 
[11] A. Dickenstein and B. Sturmfels. Elimination Theory in Codimension Two. Journal 

of Symbolic Computation (2002) 34, 119-135. 
[12] A. Dickenstein, E. M. Feichtner and B. Sturmfels. Tropical Discriminants. 
|math.AG/0510126| 

[13] D. Eisenbud and B. Sturmfels, Binomial ideals. Duke Math. J. 84 (1996). no. 1, 
1-45. 

[14] W. Fulton. Intersection Theory. Springer- Verlag, New York, 1984. 

[15] I. M. Gel'fand, M. Kapranov and A. V. Zelevinsky . Discriminants, Resultants and 
Multidimensional Determinants. Birkhauser, Boston-Basel-Berlin, 1994. 

[16] G, M. Greuel and G. Pfister. A Singular Introduction to Commutative Algebra. 
Springer- Verlag, Berlin Heidelberg New York, 2002. 

[17] I. M. Gel'fand, M. M. Kapranov, and A. V. Zelevinsky. Hypergeometric functions 
and toric varieties. Funktsional. Anal, i Prilozhen. 23(2) (1989), 12-26. 

[18] A. K. Lenstra, H. W. Lenstra, Jr., L. Lovasz, Factoring polynomials with rational 
coefficients. Math. Ann. 261 (1982), no. 4, 515-534. 

[19] H. Matsumura. Commutative Ring Theory. Cambridge studies in advanced math- 
ematics, 8. Cambridge University Press, 1989. 

[20] M. Kapranov. A characterization of A-discriminantal hypersurfaces in terms of the 
logarithmic Gauss map. Mathematische Annalen, 290, 1991, 277-285. 

[21] T. Mora, M. E. Rossi. An algorithm for the Hilbert function of a primary ideal. 
Comm. Alg. 23 (1995), 1899-1911. 

[22] M. Passare and A. Tsikh. Algebraic equations and hypergeometric series, in: The 
legacy of Niels Henrik Abel, pp. 653-672, Springer, 2004. 

[23] G.-M. Greuel, G. Pfister, and H. Schonemann. Singular 3.0, A Computer Al- 
gebra System for Polynomial Computations. Centre for Computer Algebra, 
University of Kaiserslautern, 2005, http://www.singular.uni-kl.de 



Departamento de Matematica, FCEN, Universidad de Buenos Aires, (1428) 
Buenos Aires, Argentina. 

E-mail address: alidick@dm.uba.ar, macueto@gmail.com 



